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Abstract. Cellular automata have been mainly studied on very regular graphs car- 
rying the vertices (like lines or grids) and under synchronous dynamics (all vertices 
update simultaneously). In this paper, we study how the asynchronism and the graph 
act upon the dynamics of the classical Minority rule. Minority has been well-studied 
for synchronous updates and is thus a reasonable choice to begin with. Yet, beyond 
its apparent simplicity, this rule yields complex behaviors when asynchronism is in- 
troduced. We investigate the transitory part as well as the asymptotic behavior of 
the dynamics under full asynchronism (also called sequential: only one random vertex 
updates at each time step) for several types of graphs. Such a comparative study is a 
first step in understanding how the asynchronous dynamics is linked to the topology 
(the graph). 

Previous analyses on the grid [1,2] have observed that Minority seems to induce fast 
stabilization. We investigate here this property on arbitrary graphs using tools such 
as energy, particles and random walks. We show that the worst case convergence time 
is, in fact, strongly dependent on the topology. In particular, we observe that the case 
of trees is non trivial. 

1 Introduction 

In this paper, we investigate the random process Minority: each vertex of a graph is char- 
acterized by a state or 1. At each time step (the time is discrete), one random subset 
of vertices is drawn. These vertices are updated and switch to the minority state in theirs 
neighborhood. 

Similar random processes appeared in the literature and concern different fields of appli- 
cations. For example, several studies focus on the emergence of cooperation in the iterated 
prisoner's dilemma game [3,4,5]: used as a simple social model, this field helped determine 
which ingredients are needed to foster the emergence of cooperative behaviors. The rock- 
paper-scissors game was used to model evolution of colonies of bacteria [6]. In physics, the 
Ising model (a kind of Majority rule) was introduced to study ferromagnetism (original 
paper is [7], see for instance Chapter 10 of [8] for a quick introduction): states represent 
the orientation of spin. Anti-ferromagnetism is studied by using a kind of Minority rule [9]. 
Recent works [10,11] use stochastic Minority to model the formation of quasi-crystals. 

Initially, our interest in studying Minority process comes from the field of cellular au- 
tomata (CA). CA can be seen both as a model of computation with massive parallelism and 
as a model for complex systems in nature. They have been studied with various fields of 



applications like parallel/distributed computing, physics, biology or social sciences. Most of 
the work regarding CA assumes that their dynamics is deterministic and synchronous (all 
vertices update simultaneously) and that the graph is very regular (usually a line or a 2D 
or 3D grid). Such assumptions can be questioned with regard to the applications and the 
real life constraints. Dynamics where those assumptions are perturbed have been far less 
studied and their analysis is very challenging. Here is a non-exhaustive list of related works 
about CA in literature, the models are stochastic CA since the perturbations are usually 
introduced as stochastic processes: 

— Perturbation of the updating rule: resilience to random errors [12,13,14,15,16,17], mean 
field analysis of general Markovian rules [18]. 

— Perturbation of the synchronism (i.e. of the updating scheme): empirical studies about 
resilience to asynchronism [19,20,21,22,23], mathematical analysis of some ID CA under 
full asynchronism (only one random vertex updates at each time step) or under a- 
asynchronism (each vertex updates independently with probability a) [24,25,26]. 

— Perturbation of the graph (the topology of cells): empirical studies [27,28,29], gene reg- 
ulatory networks [30,31]. 

Previous analyses focus on the effects of asynchronism on 2D Minority with Von Neu- 
mann neighborhood [1] and Moore neighborhood [2]. In this paper, we choose to investigate 
how the graph acts upon the dynamics under asynchronous updates. (Note that this is a 
special case of Interacting Particle Systems [32].) We focus on Stochastic Minority where 
the Minority rule applies to two possible states {0, 1} and under full asynchronism (at each 
time step, only one random vertex is updated with the uniform distribution). This simple 
rule already exhibits a surprisingly rich behavior as observed in [1,2] where it is studied 
for vertices assembled into a torus. Such behaviors may appear because Minority is a CA 
with negative feedback. The evolution of CA with positive feedback can be described by a 
bounded decreasing function over time [33]. Thus, the difficulty of analyzing the Minority 
rule (negative feedback) must not be confounded with the difficulty of analyzing the Ma- 
jority rule (positive feedback). Some related stochastic models like Ising models or Hopfield 
nets have been studied under asynchronous dynamics (e.g. our model of asynchronism cor- 
responds to the limit when temperature goes to in the Ising model). These models are 
acknowledged to be harder to analyze when it comes to arbitrary graphs [34,9] or negative 
feedbacks [35]. 

Let us stress that we mainly study stochastic Minority on several kinds of graphs: trees, 
cliques, bipartite graphs and compare them. General results are not precise enough to de- 
scribe its behavior in the present study. One of the aim of this paper is to show that Minor- 
ity's behavior highly depends on the topology of the graph. Our paper focus on particular 
classes of graphs. It is a first step and the reasonings may prove to be helpful to study future 
applications of Minority (as in [10,11]) and our results complete some previous results about 
Minority [1,2]. 

Here is a list of previous claims about Stochastic Minority in the literature, as well as 
some insights provided by the present paper. 

Short introduction of stochastic Minority behaviors: Figure 1 shows Minority on a 
2D grid under three different dynamics: 



— the asynchronous dynamics where each cell has a probability a to be updated indepen- 
dently from the others at each time step; 

— the synchronous dynamics where all cells are updated at each time step (a-asynchronous 
dynamics for a = 1); 

— the fully asynchronous dynamics where only one cell, randomly and uniformly chosen, 
is updated at each time step (it can be regarded as, and often behaves as the limit of 
the a-asynchronous dynamics when a tend to [25,36]). 




Fig. 1. Stochastic Minority under different a- asynchronous dynamics with N$o = 50 x 50 cells 
arranged in a 2D grid with periodic boundary condition (i.e. torus). The last column gives, for 
a £ [0,1], the empirical probability that an initial random configuration converges to a stable 
configuration before time step t s • N$o where t s = 1000 and t s = 2000 for von Neumann and Moore 
neighborhood respectively. The fully asynchronous dynamics is abusively designed by a = 0. 

Depending on the value a, Minority can exhibit two different behaviors. Experimentally, 
these phenomenon can be observed as a phase transition depending on a on the convergence 
time (see figure 1): 

— When a is almost 1 , there is a big homogeneous flashing background with some random 
noise. By flashing, we mean that all the cells of the background are not in their minority 
state and since a is almost 1, they switch their states at almost each time step. The few 
cells which are not updated create some random noise. Experimentally, the dynamics 
never reaches a stable configuration, but a highly improbable series of updated may lead 
to a stable configuration. 

— When a is almost 0, different regions made of checkerboard patterns (von Neumann 
neighborhood) or stripes (Moore neighborhood) quickly appear. The only cells which 
may switch their state are along the borders between these regions. These borders evolve 
and they eventually stabilize. Experimentally, the dynamics reaches a stable configura- 
tion in polynomial time according to the size of the configuration. 



No results were previously known for the a-asynchronous dynamics. In section 4.3, we study 
this phase transition in ID and show that it is linked to directed percolation. The proof of 
this result (theorem 34) is short because most of the technical aspects were previously done 
in [37] on a different automaton. 
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Fig. 2. Minority under fully asynchronous dynamics on 2D grids with periodic boundary condition, 
von Neumann and Moore neighborhood [1,2]. 

The previous studies [1,2] focuses on the fully asynchronous dynamics which is similar 
to the a-asynchronous dynamics when a is almost (see Figure 2). These papers analyse 
the start and the end of a classical execution on stochastic Minority on a 2D grid with von 
Neumann and Moore neighborhood. Indeed, they had to study separately the formation of 
the regions (beginning) and the evolution of the borders (ending). 

Local interactions: under fully asynchronous dynamics, patterns, which depend on the 
topology of the graph, quickly appear. In theorems 5 and 7 of [1] and theorem 9 of [2], the 
authors show that this phenomenon is due to local interactions and occurs in polynomial 
time according to the size of the configuration. When these papers were written, the authors 
did not need to be more precise but they conjectured that these patterns appear in linear 
time. Here by studying cliques (where long range interactions do not exist), we show that 
the dynamics stabilizes in linear time. This result supports the previous conjecture. 
Long range interactions: the evolution of borders between different regions often implies 
interactions between cells which could be arbitrarily far away in the graph. This long range 
interactions are harder to analyze and major differences appear between the von Neumann 
and Moore neighborhood on the 2D grid. This remark leads us to analyze different topologies. 
Here we solve a conjecture made in [1]. We exhibit biased trees were Minority behaves as a 
biased random walk and need an exponential time to converge toward a stable configuration. 
In the previously considered topologies, Minority always converges in polynomial time under 
fully asynchronous dynamics. 

Capacity of simulation: one important aspect of Minority is the diversity of phenomena 
embedded in this rule. On a-asynchronous dynamics, there is a phase transition between 
two different behaviors. In section 4.3, we establish a link between this phase transition and 
directed percolation. Previous studies [1,2] have shown that the fully asynchronous dynamics 
occurs in two steps. We show here (Section 4.1) that the first step acts as a coupon collector 
on cliques. When long range interactions occur, the dynamics may behave as competition 



between 2D regions, different kinds of random walks. Other behaviors may be encoded on 
stochastic Minority. It is surprising that a simple rule may produce such different behaviors 
from random configurations on regular topologies. 

Bipartite graphs: one aim of this paper is to generalize tools previously designed to study 
Minority on grid. Basic tools may be generalized to any graph but we realized that advanced 
tools can only be generalized to bipartite graphs. This leads to an explanation between the 
difference of behaviors previously observed on 2D grids. Note that even if our tools are 
helpful to analyze bipartite graphs, various and complicated behaviors may already appears 
in this class of graphs. 

2 Model 

In this paper we consider Stochastic Minority on arbitrary undirected graphs. 

Definition 1 (Configuration). Let Q = (V, £) be a finite undirected graph with vertices V 
and edges 8 . Q = {0, 1} is the set of states (0 stands for white and 1 stands for black). The 
vertices are also called cells and N := \V\ denotes their number. The neighborhood Mi of 
a vertex i is the set of its adjacent vertices (including itself). A configuration is a function 
c : V — >• Q fci denotes the state of vertex i in configuration c). 

Definition 2 (Stochastic Minority). We consider the following dynamics 5 that asso- 
ciates with each configuration c a random configuration c' obtained as follows: a vertex i G V 
is selected uniformly at random ( we say that vertex i is fired ) and its state is updated to the 
minority state among its neighborhood (no change in case of equality), while all the other 
vertices remain in their current state. Formally: 



and (6(c)) k = cj~ for all k ^ i. In a configuration, a vertex is said to be active if its state 
changes when the vertex is fired. The random variable d denotes the configuration obtained 
from an initial configuration c, after t steps of the dynamics: c° = c and d = 5 t (c) for all 
t ^ 1. (The notation 5 t (c) means 5 o ... o 5(c), t times). 

Definition 3 (Attractors and limit set). For the dynamics induced by 5, a set of con- 
figurations A is an attractor if for all c, d G A, the time to reach d starting from c is finite 
almost surely. In the transition graph where vertices are all the possible configurations and 
arcs (c, d) satisfy F (5(c) = d) > 0, attractors are the strongly connected components with 
no arc leaving the component. The union of all attractors is denoted A and called the limit 
set. 

Definition 4 (Convergence and hitting time). We say that the dynamics 5 converges 
from an initial configuration c° to an attractor A (resp. the limit set A) if the random 
variable T = min{t | d G A} (resp. T = mm{t | d G A}) is almost surely finite. 
The variable T is a hitting time. 
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Since we only consider finite graphs, the dynamics 5 always converges from any initial 
configuration to A, and T is well-defined (an exception is section 4.3, where we consider an 
infinite graph and discuss the convergence to an attractor). 

A hitting time T is defined for a given graph and a given initial configuration. We are 
interested in the worst case (i.e. largest hitting time) among all possible initial configurations 
and among all graphs of a given class of graphs. 

3 Tools 

3.1 Energy, Potential 

As in the Ising model [9] or in Hopfield networks [35], we define a natural global parameter 
similar to an energy: it counts the number of interactions between neighboring vertices in 
the same state. This parameter will provide key insights into the evolution of the system. 

Definition 5 (Potential). The potential V{ of vertex i is the number of its neighbors (in- 
cluding itself) in the same state as itself. If V{ ^ then the vertex is in the minority 
state and is thus inactive; whereas, if V{ > then the vertex is active. A configuration c 
is stable if and only if for all vertex i G V , V{ ^ 

Definition 6 (Energy). The energy of configuration c is E(c) = J2iev( v i ~ !)• 

The energy of a configuration is always non-negative. There are configurations of energy 
if and only if Q is bipartite: those stable configurations are the 2-colorings of Q. More 
generally, we have: 

Proposition 7 (Energy bounds). The energy E satisfies 2\£\ — 2C max ^ E ^ 2\£\, where 
C max is the maximum number of edges in a cut of Q . 

Proof. The bounds are direct consequences of the definitions: for any configuration, E = 
2\£\ - 2\C\ where C is the cut {{i, j} G £ \ c { = and cj = 1}. 

As a consequence, computing the minimum energy for arbitrary graphs is NP-hard: it is 
equivalent to computing MAXCUT. 

Lemma 8 (Energy is non-increasing). The energy is a non- increasing function of time 
and decreases each time a vertex i with potential stricly larger than fires. 

Proof. When an active vertex i of potential V{ fires, its potential becomes |A/";| — vi + 1, and 
the energy of the configuration becomes E + 2\Mi\ — 4^ + 2. □ 

Corollary 9. A configuration c belongs to the limit set if and only if no sequence of updates 
would lead the energy to decrease, i.e. if and only ifVt, P (E^) > E(S(c t )) | c° = c) =0. 

Proof. If the energy decreases when updating c to c ; , then c will never be reached again 
(because energy is non- increasing) . Reciprocally, any update that keeps the energy constant 
is reversible: the fired vertex can be fired again to get back to the previous configuration. □ 



Remark 10. Since firing a vertex of odd degree makes the energy decrease, such vertices are 
inactive in the limit set. 

Definition 11 (Particle). Let c be a configuration on Q — (V,£), an edge {i,j} holds a 
particle if Ci = Cj. A configuration is fully characterized (up to the black/white symmetry) 
by its set of particles located atV^S. 

(Note that the converse proposition "any subset V C £ corresponds to a configuration" is 
true if and only if the graph is a tree). 

The energy of a configuration is clearly equal to twice its number of particles. With the 
particle point of view, when firing a vertex i of degree deg(i) = \J\fi\ — 1, if the number of 
incident edges holding a particle is a least deg ^ ? these particles disappear but new parti- 
cles appear on the incident edges (if any) which had no particle (as illustrated on Fig. 3). 
Otherwise the particles do not move. 




Fig. 3. Transfers of particles (diamonds) when firing vertex i. 



Switching between the coloring and the particle points of view may simplify the descrip- 
tion of the configurations and the dynamics, e.g. when the energy is low and the dynamics 
comes to random walks of a few particles. 

3.2 Bipartite Graphs 

A graph is bipartite when its vertices can be partitioned into two sets such that every edge 
goes from one set to the other (or equivalently, when it is 2-colorable) . Bipartite graphs allows 
us to use another, easier point of view (the dual configuration) and to easily determine if a 
configuration is in the limit set. 

Dual configuration We now introduce dual configurations as in [1] (section 3.2), and their 
dual rule to facilitate the study of the dynamics on trees. In this dual dynamics equivalent to 
Stochastic Minority, the stable configurations of minimum energy are the two configurations 
all black and all white, and the regions which compete are all white versus all black subtrees. 

The dual rule is almost a majority rule, but in case of equality among the neighbors of a 
vertex, the state of this vertex is flipped each time it is updated. This "instability" prevents 
many results about majority rules to apply to our case. 

Definition 12 (Dual configurations). Consider a graph Q and fix a vertex r (the "root"). 
For any configuration c on Q , its dual configuration c is defined as Ci = q if hi is even and 
&i = 1 — Ci if hi is odd, where hi is the distance from r to i (see Fig. 4)- The mapping c ^ c 
is a bisection on the set of all the configurations; more precisely c = c. 



An equivalent definition consists in making a XOR with the 2-coloring of Q such that r 
is white. The duals of the configurations of minimum energy are the configuration where 
all vertices are black or all vertices are white. 



Q p 




Fig. 4. A configuration c and its dual configuration c (with regard to root r). 



Proposition 13 (Dual dynamics). Consider a sequence (c l ) for the Stochastic Minority 
dynamics 5 and the sequence (c?) of the dual configurations, and define the dual dynamics 
5 as 5(c) = 6(c) so that c £+1 = £(c £ ). Then the dynamics 6 is also a stochastic CA. It asso- 
ciates with each configuration c a random configuration c' by updating one random vertex i 
uniformly with the rule which selects the majority state in the neighborhood of i excluding 
itself (in case of equality its state changes): 

g/ = f 1 if E^Mxw £ j > ^^2^ or (E^Mxw £ j = and ^ = °) 

Ci I if E^MxW * < ^ or E Ms{i) c 3 = J^=i and c % = 1) 

By construction, the dual sequences (c £ ) and (c £ ) as well as their corresponding dynamics S 
and S are stochastically coupled (see [38]) by firing the same random vertex at each time 
step. 

Definition 14 (Dual potential & Energy). The dual potential V{ of vertex i is the 
number of its neighbors (excluding itself) in a different state than itself Iffy < ^j -1 then 
the vertex is in the majority state and is thus inactive; whereas, if fy ^ ^j" 1 ^ en ^ e 
vertex is active. The dual energy E is the sum of the dual potentials over all the vertices. 

Given a configuration c and its dual c, the potential of any vertex i in c is equal to the 
dual potential of vertex i in c plus 1. Thus the dual energy of c is exactly the energy of c. 

Theorem 15. Consider a configuration c and its dual c. Then, if the same vertex fires, 
5(c) = (5(c). We thus have the following commutative diagram: 



original c ^ c ' 



Proof. Consider a configuration c and its dual c. A vertex q of c is active if and only if 
v i ^ Since ^ = # + 1 and ^ is an integer, a vertex C{ of c is active if and only if 
Vi < • A ^ 2 ~ 1 : that is to say vertex q of c is active. According to definition 12, if the same 
vertex fires and a vertex of c is active if and only if the corresponding vertex of c is active 
then 5(c) = 5(c). □ 

Distance to a Stable Configuration In this section we describe an algorithm (algo- 
rithm 1) that gives a sequence of updates that leads to the limit set. It is then easy to test 
for the limit set: the input configuration belongs to the limit set if and only if the energy is 
the same between the input and output configurations. 



Fact 16 An attractor A decomposes the graph into three sets of vertices: 

1. the vertices that are in the state for every configuration of A; 

2. the vertices that are in the state 1 for every configuration of A; 

3. the vertices that can be either in the state or 1, depending on the configuration in A. 



Algorithm 1: Membership to the limit set: check that E(c') = E(c). 
Input: A configuration c. 

1 while There is an active white vertex i do Fire i (i becomes black) 

2 while There is an active black vertex i do Fire i 

3 while There is an active white vertex i do Fire i 
Output: The configuration c' at the end of phase 2. 



Proposition 17. The configuration d returned by algorithm 1 is in the limit set. 

Corollary 18. The input configuration c is in the limit set if and only if the energy has not 
decreased during execution of the algorithm. 

Proof. We first prove that the vertices in state (white) at the end of phase 1 of the 
algorithm cannot switch to state 1, whatever the sequence of updates, i.e. they are in Case 1 
of Fact 16. Indeed, assume instead that there exists a vertex i and sequence of configurations 
c 1 , c 2 , c k such that 

— c 1 is the configuration at the end of phase 1; 

— each configuration is the result of firing one vertex in the previous configuration; 

— cj = and c\ = 1. 

Let c £ be the configuration just before the first update of this sequence that fires an active 
vertex j in state (there exists one since at least i will be fired in this sequence). But j 
must have been already active in c 1 , since it had at least as many neighbors in state 1 as 
in c £ (this is a monotonicity argument). The algorithm thus would not have exited the first 
"while" loop, which is a contradiction. 

Same holds for the vertices in state at the end of phase 3. Also, by symmetry, the 
vertices in state 1 at the end of phase 2 are in Case 2 of Fact 16. 



Finally, since the remaining vertices were white at the end of phase 2, they can be made 
white by a sequence of updates starting from the configuration c 1 (the configuration at the 
end of phase 1). By monotonicity, they can be made white by the same sequence of updates, 
starting from any configuration reachable from c 1 . This means that the configuration where 
all those remaining vertices are white is always reachable, and is thus in the attract or. This 
configuration is the one at the end of phase 2. Which concludes the proof. 

By symmetry, those remaining vertices can be made black from the configuration at the 
end of phase 2, which is in the attractor, so they are in Case 3 of Fact 16. □ 

Termination of Algorithm 1 is clear: phase 1 increases the number of vertices in state 1 
at each iteration and is thus completed in O(N) iterations. The same remark applies to the 
other phases. 

Three phases are necessary to correctly classify the vertices into the three cases of fact 16, 
as shown on Fig. 5. Two phases are enough to reach the limit set. 




Starting from this configuration, phase 1 of the algorithm makes 
the vertex a become black. Then, phase 2 makes its left neighbor 
b become white, as well as a. 

b is now definitely white (Case 1 of Fact 16), but was not iden- 
tified as such at the end of phase 1. Thus, a third phase is 
necessary with this algorithm to identify b. 

No other vertex will be active during the execution of the algo- 
rithm. 



Fig. 5. Three phases are necessary. 



Stable Configurations Stable configurations on trees are characterized in the appendix, 
Section A. 

3.3 Non Bipartite Graphs 

It seems that the only in-depth study of a non-bipartite graph is the analysis of the 2D grid 
with Moore neighborhood and periodic boundary condition [2]. This study is harder and 
more complicated than any bipartite graphs considered so far. This complexity appears in 
the structure of the stable configurations. For example, Figure 6 shows a stable configura- 
tion, illustrating many ways in which the four different stable patterns can be intricated. 
Perturbing such stable configurations leads to original random walks (see Figure 7). Parti- 
cles appear and move along the borders. They deform the border when they move. When 
two particles collide, they disappear. When borders are too perturbed, the whole structure 
of the stable configuration collapses. As opposed to the bipartite case, we are currently not 
able to analyze such dynamics and compute for any configuration a sequence of updates 
leading to a stable configuration. 



Fig. 6. A complicated 16 x 32 stable configuration. There are four stable patterns (V,V',H,H') and 
different types of junction between the different regions(A, A', B, B' and C). 




7. a - A sequence of updates in a configuration 7.b - A sequence of updates where the rails 
starting with 4 particles where two of them move cannot sustain the perturbations due to the 
along rails and ultimately vanish after colliding movements of the particles: at some point, 
with each other. rails get to close with each other, new active 

cells appear, and part of the rail network col- 
lapses. 

Fig. 7. Some examples of the complex behavior of particles in a 20 x 20 configuration. 



4 Behaviors 

4.1 Decreasing Energy (Clique, cycle and paths) 

We prove in this part that Stochastic Minority on cliques behaves as a coupon 
collector (see [39], page 210). This easy result implies a fast convergence to 
the limit set. 

Theorem 19. Stochastic Minority on cliques hits the limit set after 0(N log N) steps on 
expectation. If N is even, the (^) attractors are the half black and half white configurations, 

each one stable. If N is odd, the only attractor is the set of 2( n^i ) configurations having 

2 

one more (resp less) black than white vertices. 




Proof. Let rib be the number of black vertices of a configuration. Since the neighborhood of 
a vertex is V, the potential of a black vertex is n^. If rib > ^±^- (resp. rib < ^r^) then firing 
a black (resp. white) vertex decreases the energy and the configuration is not in A. Thus a 
configuration in A must have ^ rib ^ Consider such a configuration: 

— If N is even then all vertices have potential y and these (^) configurations are stable. 

— If N is odd we call C5 (resp. C w ) the set of configurations where rib = (resp. rib = 
^Y^). White (resp. black) vertices of a configuration in C5 (resp. C w ) are inactive and 
black (resp. white) vertices are active, firing one of them leads with constant energy to 
a configuration of C w (resp. Cb). Thus from any configuration of C w U there is no 
sequence of updates that causes a drop of energy and C w U C5 = A. 

Now, we prove that A is made of only one attractor. Let d(c,c f ) := | {i | q ^ c • } | be 
the distance between two configurations. Consider in A two configurations c ^ c' . By 
symmetry we can assume c G C5. Since c has at least as many black vertices as c', there 
is a vertex i black in c and white in d . Firing i decreases the distance. Iterating this 
argument, one finds a path from c to d in A. 

Now consider a configuration where rib > ^2^- As long as the configuration does not 
belong to A, the white vertices are inactive. When |"n& — ^jr^] black vertices have fired, the 
configuration is in A. At each time step there is a probability ^ to fire a black vertex. This 
kind of dynamics is known as coupon collector and T = 0(N\ogN). □ 

Moreover, it is easy to see that in the odd TV case, the attractor has a simple structure: 
this is a bipartite graph, composed of configurations with white vertices on one side, 
and configurations with white vertices on the other. 

Cycles and paths Cycles and paths forms the class of connected graphs of maximum 
degree 2. 

On cycles and paths, the particle point of view is convenient and one can prove 
that Stochastic Minority behaves as random walks of annihilating particles on 
a discrete ring (see [24,40]). On the right, the particles are the diamonds. 

Theorem 20. Stochastic Minority on cycles and paths hits the limit set after 0(N 3 ) steps 
on expectation. If N is even, the two attractors are the 2-colorings of the cycle. If N is odd, 
the single attractor is a cycle in the transition graph composed of all the configurations with 
only one particle. 

Proof. The movement of particles provide a nice framework for the analysis. Firing a vertex 
incident to a particle either attract the particle on the next edge if it is free, or annihilate 
both particles on each side of the vertex. Consequently, the dynamics boils down to the 
analysis of random walks of annihilating particles on a discrete ring. If the number N of 
vertices is even (resp. odd), any configuration has necessarily an even (resp. odd) number of 
particles. This number decreases by annihilations until there is no (resp. only one) particle. 




The attract or is reached at this point since the energy cannot decrease further. This proof 
can be easily adapted to graphs which are paths. 

To bound the expected hitting time of the limit set, associate with each configuration c l 
a weight X t which is the maximum distance between two consecutive particles if there are 
at least two particles, or N if there is only one particle, or N + 1 if there is no particle. For 
all t, X t e {1, • • • , N + 1} and c t belongs to the limit set if and only if X t e {N, N + 1}. Let 
AX t+ i = X t +i — X t . One can check that ~E(AX t +i | c t = c) ^ for any configuration c. 
Moreover E ((/AX^+i) 2 | c t = c) ^ 3/N for any configuration c not in the limit set. Conse- 
quently X} - 3t/N is a sub-martingale and we can apply the Stopping Time Theorem to the 
stopping time T = min {t ^ | X t e {N, N + 1} }. It gives E (T) = 0(N 3 ) which is thus an 
upper bound on the expected hitting time of the limit set. 

This proof applies with no modification to the hitting time on graphs which are paths. 

□ 



4.2 Long Range Interaction and Exponential Convergence (Trees) 

In this part, we introduce biased trees (Definition 24 and Figure 9) such that the dynamics 
c converges in exponential time on this topology (Theorem 31). Vertices of biased trees 
have degree at most 4. In fact, biased trees simulate biased random walks (Definition 21) 
which converge in exponential time. Biased trees are created from small trees called widgets 
(Definition 23 and Figure 8) arranged on a line. Except at the ends, this line of widgets is 
made of "gates" . According to the configuration, these gates are either locked, unlocked or 
stable (Definition 25). On a correct configuration (Definition 26), the line of gates is split into 
two regions: all gates on the left side are stable and all gates on the right side are unstable 
(locked or unlocked). In a correct configuration, three different events may be triggered with 
the same probability 1/N (Fact 27 and Corollary 30): 

— the rightmost stable gate becomes an unlocked gate; 

— the leftmost unstable gate becomes stable if it is unlocked; 

— the leftmost unstable gate is switched from locked to unlocked or the contrary. 

Thus stable gates tend to disappear. This dynamics will ultimately converge to the stable 
configuration df (Definition 28). To reach this configuration all gates must be stable. Thus 
it takes an exponential time for the dynamics c to converge on a biased tree with an initial 
correct configuration. 

Definition 21 (Biased random walks). A Biased Random Walk is a sequence of random 
variables (J*Q)^ defined on {0, ... , n} such that for all i ^ 0: 

— P (X i+1 = l\Xi=0) = l (reflecting barrier at 0). 

— P (Xi+i = n | Xi = n) = 1 (absorbing barrier at n). 

) (X i+1 =x-l\X i = x) 

-3a,6eM + \fx e {l,...,n- 1} { +P (X i+1 = x + 1 1 X { = x) = 1 

< a < P(X i+ i = x + 1 \Xi = x) < b < 1/2 



Theorem 22. Let T := min{z ^ | Xi = n} be the absorption time at n and for all ^ 

k ^ n, let Efc (T) := E (T | Xq = k) be its expectation starting from k. Then 



9 k (b) < E fe (T) < fc (o) 
where k (p) = ((^)" - - £ 



n —k 
-2p' 



This theorem is a direct consequence of classical analysis of random walks on {0, 



i} where 



- F(X i + 1 = l\X i = 0) = l 

- ¥(X i+1 = n\Xi = n) = 1 

- Vxe {!,..., n-l} /P(^i=- + l|^=»)=P 
Solving the following system of equations [39]: 



with p + q = 1 . 



E k (T) = p(l + E fc+1 (T)) + (1 - p)(l + E fe _i (T)) Vfc e {0, n} 

E n (T) = 

E (T) = l + Ei (T) 



one 



gets E fc (T) = M^-P) _ _ 



Definition 23 (Widgets). ^4 Widget W is a tree T = (V, £,&) where b e V is called the 
bridge. We consider the three widgets described in Figure 8. a: head, gate and tail and the 
three configurations q,c u ,c s for gates. 



Definition 24 (Biased trees). Let (Wi)o^i^ n +i be a finite sequence of widgets where 
W{ = (Vi,£i,bi). From this sequence, we define the tree T = (V, £) where V = U^j^H and 
£ = (^7=0^) U(^?=o^^+i)- Abusively we also denote by (Wi)o^i^n+i ^ £ree generated by 
this sequence. A biased tree of size n is a finite sequence of widgets (Wi)o^i^n+i where Wo 
is a head, for 1 ^ i ^ n, W{ is a gate and W n +i is a tail. 



Definition 25 (Stable and unstable gates). Consider a biased tree (W^o^i^n+i an d a 
configuration c. We denote by cw if the restriction of c to widget W{. We say that gate i is 
locked if = Q, unlocked if = c u and stable if = c 8 . An unstable gate is a gate 
which is locked or unlocked. 

Definition 26 (Correct configuration). Configuration c is correct if vertices of the head 
are black, the tail is white, and there exists a j such that for all 1 ^ i ^ j gate i is stable 
and for all j < k ^ n gate k is unstable. We say that configuration c is on position j. We 
denote by Pos(c) the position of configuration c. The position is unlocked if j = n or gate 
j + 1 is unlocked, the position is locked otherwise. 




Unstable Gates 

8. a - The 3 widgets used in the con- 8.b - The three configurations ci,c u and c s . 

struct ion of a biased tree. Gray denotes 
the fact that the vertex state is not rep- 
resented. 

Fig. 8. Widgets used in the construction of a biased tree. 

Fact 27 (Active vertices) Consider a correct configuration c on position j. The active 
vertices of c are: 

— Vertex bj if j ^ 0. 

— Vertex bj + \ if j ^ n and gate Wj+\ is unlocked. 

— Vertex fi if j < i ^ n. 

— Vertex b n +i if j = n. 

Proof. Consider a correct configuration c on position j. The only vertices which may be 
active in c are vertices bi and fi for 1 ^ i ^ n and vertex 6 n +i. Vertex 6 n+ i is active if and 
only if c(b n ) = 1 that is to say gate W n is stable. For all 1 ^ i ^ n, vertex fi is active if 
and only if the gate W{ is unstable that is to say j < i $J n. For all 1 ^ i ^ n, vertex bi 
is inactive if c(^_i) = c(bi) = c(^ + i). Thus among vertices (6i)i^^ n , only vertices bj and 
bj+i may be active: vertex bj is active and vertex bj+\ is active if gate Wj+i is unlocked. □ 

Definition 28 (Final configuration). The final configuration Cf is the configurations 
where vertices of the head are black, the tail is black and every gate is stable. We say that 
Cf is on position n + 1, Pos(df) = n + 1. 

Lemma 29. Configuration Cf is stable. 

Proof. Consider the correct configuration c on position n. According to fact 27, only vertices 
b n and 6 n+ i are active. If vertex 6 n+i fires, these two vertices become inactive and no other 
vertex becomes active. Firing vertex 6 n+ i leads to configuration Cf. Thus Cf is stable. □ 

Corollary 30. Consider a correct configuration c, then configuration c' = 5(c) is either 
correct or Cf. Moreover \Pos(c') — Pos(c)\ ^ 1. 




O : active cell 

Fig. 9. A biased tree and a correct configuration on unlocked position 2. 



Proof. Consider a correct configuration c on position j and the configuration c! = 5(c). If 
an inactive vertex fires then c' = c. Now consider that an active vertex fires (see fact 27): 

— if j and vertex bj fires: then gate Wj becomes unlocked and c! is a correct configu- 
ration on unlocked position j — 1. 

— if j 7^ n and vertex fires: then gate Wj+i becomes stable and c! is a correct 
configuration on position j + 1. 

— if vertex fi fires with j < i ^ n: then gate Wi becomes unlocked (resp. locked) in cf if it 
is locked (resp. unlocked) in c. Configuration c' stays correct and on position j. 

— if j = n and vertex b n+ i fires: then c! = df. □ 

Theorem 31. On biased trees of size n (i.e. N = 8n + 4 vertices), starting from a correct 
configuration, Stochastic Minority converges almost surely to c? . Moreover the hitting time T 
of the limit set satisfies <9(1.5 n ) ^ E (T) < <9(n4 n ). 

Proof. Consider a biased tree of size n, an initial correct configuration c° on position 
and the sequence (c t ) t ^o. Dynamics 5 converges almost surely from initial configura- 
tion c° and c T = Cf. We define the sequence of random variable (U)i^o as to = and 
ti+\ = min {t > U | Pos(c ti+1 ) ^ Pos(c ti ) or Pos(c ti+1 ) = n + 1 }. Consider the sequence of 
random variable (J*Q)^o such that = Pos(c ti ). According to corollary 30, Xi\ = 1. 

Consider a configuration c* on locked position n > j > then firing vertex 6j leads to 
a configuration on position j — 1 and firing vertex fj leads to a configuration on unlocked 
position j. Firing other vertices does not affect the position of the configuration. Consider a 
configuration & on unlocked position n > j > then firing vertex bj leads to a configuration 
on position j — 1, firing vertex fj leads to a configuration on locked position j and firing 
vertex bj+\ leads to a configuration on position j + 1. Firing other vertices does not affect 
the position of the configuration. A vertex has a probability 1/N to fire where N = 4 + 8n. 
Thus, the evolution of a configuration on position < j < n can be summarized as: 



(N-2)/N 



A basic analysis yields that: 

- if 1 < x < n then P (X i+1 = x + l\ Xi = x) = 1 - W [X i+1 = x - 1 | X { = x) and 1/5 < 
P(xi+i = x + 1 |X f = x) < 2/5. 

- P(X i+ i = l|X i = 0) = l. 

- P(X i+ i =n+l|Xi = n+l) = 1. 

Thus the behavior of (J*Q)^ is as described in definition 21. We define the random vari- 
able T" = min {i ^ |X^=n} which corresponds to the first time when all gates are stable, 
then 6>(§ n ) < E(T') < <9(4 n ) (see Def. 21). We call c /_1 the correct configuration on posi- 
tion n (i.e. all gates are stable). Then c T = c*^, c T_1 = c*^ -1 and P (c £+1 = c*^ | c £ = c*^ -1 ) = 

1/2. Thus, E (T) = <9(E(V))- By definition, i r , = = E£i[(*i-*i-i)l**<r']- 

Since there are at most 2 vertices which may modify the position of a correct configuration, 
we have 1 < E(t i+1 - < 0(n). Thus E"i(1^<t0 < E (t T >) < e(£ t ~i(ralt 4 <T'))- We 
conclude that 6>((§) n ) < E (T) < 6>(n4 n ). □ 

Subcase: Binary Trees Converge Rapidly In this section, we note that on binary 
trees, i.e. trees where the degrees are at most 3, the dynamics ends by fixing the states of 
the vertices of degree 1 and 3 (see Remark 10) and some isolated particles may remain and 
oscillate on disjoint paths. 

Definition 32 (Path). 

In this subcase, we call path a connected subgraph where all nodes Q 
but the end nodes have degree 2 (end nodes must thus have degree 1 • — • — • — • — • 
or 3). An example is composed of the black nodes on this figure: ^ 

Theorem 33. Stochastic Minority on trees with degrees at most 3 hits the limit set in 0(N A ) 
steps on expectation. The attractors of a tree T are in bisection with the matchings of the 
reduced tree T' where each path of T has been replaced by an edge, then each leaf has been 
removed. 

Proof. We study here the movements of particles in the initial tree T. One can divide T 
into its induced subgraphs which are paths. Those paths link the vertices of odd degree (1 
or 3). The reduced tree T' is obtained by replacing each path by an edge, then removing 
the leaves. 

Consider a configuration c on T which belongs to an attractor. There cannot be two 
particles on the same path , otherwise a sequence of updates could lead to the collision of 



these two particles and thus to an energy decrease. In the same way, there cannot be two 
particles on two paths which share a common extremity i. This extremity would necessarily 
be a vertex of degree 3, then a sequence of updates could position the two particles on the 
edges incident to i. Firing i at that time would decrease the number of particles by at least 1 
and lead to an energy decrease. Finally, there cannot be a particle on a path which has an 
end of degree one, since the particle could disappear at this end. 

Reciprocally, for configurations where there are no particles on 
the same path nor on adjacent paths nor on a path with an end of 
degree 1, the number of particles on each path is constant. This 
proves that the energy cannot decrease, and that such configura- 
tions belong to the limit set. This also establishes the bijection 
between attractors and matchings of T' (the matching indicates 
where the isolated particles are located in T). An illustration of 
T' is on the opposite figure, particles are next to the edges. 

To prove the bound on the expected hitting time of the limit 
set, we find a bound on the time until at least one particle disap- 
pears. Consider a configuration where there exist two particles on a same path of length n. 
One can suppose that these two particles follow a random walk on this path with reflecting 
barriers at each extremity, unless they collide with another particle (leading to the loss of 
two particles in the tree) or unless one of the two particles leaves the path (leading to the 
loss of at least one particle in the tree since leaving necessarily involves a vertex of degree 3 
fired with two incident particles). Thus under the condition that they have not disappeared 
before, a bound on the expected time elapsed until they collide can be derived from classical 
studies of random walks with reflecting barriers [39]: this expected time is bounded by 0(n 3 ), 
and consequently it is also a bound on the time until at least one particle disappear. 

Consider a configuration where there exist two particles on 
two paths of respective lengths n and m, sharing a common 
extremity: the vertex i. An illustration is on the opposite figure. 
With the same reasoning, assuming that the two particles have 
not led to the removal of another particle means that they follow 
random walks on their respective paths with reflecting barriers 
at the extremities. Then they can only disappear by being both 
incident to vertex i when vertex i is fired. By analyzing the 
two-dimensional random walk corresponding to the evolution of 
the respective distances to vertex z, this event occurs after at 
most 0(max(n, m) 3 ) steps on expectation, as proved in Section B 
(or in [41]) using standard tools for multi-dimensional random walks. 

Finally, from any configuration which does not belong to the limit set, at least one particle 
disappear within 0(N 3 ) steps on expectation. Since the number of particles in any initial 
configuration on a tree is bounded by AT, the expected time to hit the limit set is bounded 
byO(A- 4 ). * □ 





4.3 Phase transition 



In this part, we consider the infinite graph where the set of vertices is N and vertex i has 
two neighbors: i — 1 and i + 1. We consider the initial configuration c %n%t where the state of 
vertex i is 1 if i = and z mod 2 otherwise. This configuration possesses only three vertices 
which are not in their minority state: —1, and 1. Updating vertex leads to a stable 
configuration. 

We consider the fully asynchronous dynamics: at each time step, only one vertex is 
updated and this vertex is selected uniformly at random among the active vertices. Note 
that the set of active vertices is always finite, so this random selection is feasible. The limit 
set of an execution of Minority starting from the initial configuration c l7llt is composed of 
the stable configuration where the state of vertex i is i mod 2. 

We denote by P a (c* m£ ) the probability that the dynamics reaches a stable configuration 
from the initial configuration c init . If ¥ a (c m * £ ) < 1 then the expected time to reach the 
limit set is infinite. Experimentally, there is a phase transition on F a (c mU ) depending on a 
with a critical value a c « 0.5 such that: 

- if a < a c then F a (c init ) = 1; 

- if a > a c then F a (c init ) < 1 

Our result link the critical value a c to the critical value 0.6298 < p c < 2/3 of directed 
percolation. 

Theorem 34. If a ^ $/l - (1 -p c ) 4 then F a (c init ) < 1. The quantity $/l - (1 -p c ) 4 is 
in [0.993; 0.996]; thus if a ^ 0.996 then F a (c init ) < 1. 

Proof In [37], this theorem is proved for this dynamics with the rule FLIP-IF-NOT-EQUAL 
(an updated vertex switches its state if at least one of its neighbor is in a different state 
than itself) from the initial configuration c' init where vertex i is in state if i ^ and 1 
otherwise. We simply prove here that flip-if-not-equal is the dual dynamics of Minority 
obtained by switching state of vertices i if i mod 2 = 1: 

— if vertex i is active for Minority then at least one if its neighbors is in the same state 
as itself. Thus in the dual configuration, vertex i has at least one of its neighbors in a 
different state as itself. 

— if vertex i is inactive for Minority then both its neighbors are in a different state than 
itself. Thus, in the dual configuration, vertex i has both its neighbors in the same state 
as itself. 

Also, c' l7llt is the dual configuration of c l7llt . 
5 Conclusion 

The table below sums up the different worst case average hitting times of a limit set for 
different topologies and compares the fully asynchronous dynamics to the synchronous one. 
In the case of the torus under fully asynchronous dynamics, it is conjectured that this average 
"convergence" time admits a polynomial bound in the number of cells. 





Fully Asynchronous 


Synchronous 


Path or Cycle 


Poly 


Exp [33] 


Tree, max degree ^ 3 


Poly 


Exp [33] 


Tree, max degree ^ 4 


Exp 


Exp [33] 


Torus, von Neumann neighborhood 


Poly ? [1] 


Exp [33] 


Torus, Moore neighborhood 


Poly ? [2] 


Exp [33] 


Clique 


Poly 


Poly [33] 



The Minority rule admits a rich range of behaviors under full asynchronism. The case 
of trees has shown that the average hitting time of limit sets is not necessarily polynomial 
under full asynchronism (there is a threshold on the maximum degree). For now, it is an 
open problem to predict from the graph topology whether the dynamics will converge fast 
or slowly. Following this work, a challenge is to identify the graph parameters that guide 
this convergence speed, as well as extending such results to other updating rules. 

Acknowledgements. We are very grateful to the referees for their comments and sugges- 
tions which have improved the paper. 

References 

1. Regnault, D., Schabanel, N., Thierry, E.: Progresses in the analysis of stochastic 2D cellular 
automata: A study of asynchronous 2D minority. Theoretical Computer Science 410 (47-49) 
(2009) 4844-4855 

2. Regnault, D., Schabanel, N., Thierry, E.: On the analysis of "simple" 2D cellular automata. Dis- 
crete Mathematics and Theoretical Computer Science (DMTCS) 12(2) (2010) 263-294 online 
publication, http:/ /www. dmtcs.org/dmtcs-ojs/index.php/dmtcs/article/view/1337. 

3. Kittock, J.E.: Emergent conventions and the structure of multi-agent systems. In: Proceedings 
of the 1993 Santa Fe Institute Complex Systems Summer School. (1993) 

4. Dyer, M., Goldberg, L.A., Greenhill, C, Istrate, C, Jerrum, M.: Convergence of the iterated 
prisoner's dilemna game. Combinatorics, Probability and Computing 11 (2002) 125-147 

5. Mossel, E., Roch, S.: Slow emergence of cooperation for win-stay lose-shift on trees. Machine 
Learning 67(1-2) (2006) 7-22 

6. Kerr, B., Riley, M.A., Feldman, M.W., Bohannan, B.J.M.: Local dispersal promotes biodiversity 
in a real-life game of rock-paper-scissors. Nature 418 (2002) 171-174 

7. Ising, E.: Beitrag zur Theorie des Ferromagnetismus. Zeitschrift fiir Physik 31 (1925) 253-258 

8. Srivastava, R., Ashok, J.: Statistichal Mechanics. Prentice-Hall of India (2005) 

9. McCoy, B., Wu, T.T.: The two-dimensional Ising model. Harvard University Press (1974) 

10. Bodini, O., Fernique, T., Regnault, D.: Quasicrystallization by stochastic flips. Proceedings of 
Aperiodic 09 (2009) to appear 

11. Bodini, O., Fernique, T., Regnault, D.: Stochastic flips on two- letter words. 
In: Proceedings of the Seventh Workshop on Analytic Algorithmics and Com- 
binatorics (ANALCO 2010), SIAM (2010) pages 48-55 , online publication 
http: / / www.siam.org/ proceedings / analco /2010/ analcolO.php. 

12. Gacs, P.: Reliable computation with cellular automata. Journal of Computer and System 
Sciences 32(1) (1986) 15-78 

13. Gacs, P.: Reliable cellular automata with self-organization. Journal of Statistical Physics 
103(1/2) (2001) 45-267 



14. Gacs, P., Reif, J.: A simple three-dimensional real-time reliable cellular array. Journal of 
Computer and System Sciences 36(2) (1988) 125-147 

15. Toom, A.: Nonergodic multidimensional systems of automata. Problems of Information Trans- 
mission 10 (1974) 239-246 

16. Toom, A.: Stable and attractive trajectories in multicomponent systems. Advances in Proba- 
bility 6 (1980) 549-575 

17. Toom, A., Vasilyev, N., Stavskaya, O., Mityushin, L., Kurdyumov, G., Pirogov, S.: Discrete 
local Markov systems. In: Stochastic cellular systems: ergodicity, memory, morphogenesis. 
Manchester University Press (1990) 1-182 

18. Balister, P., Bollobas, B., Kozma, R.: Large deviations for mean fields models of probabilistic 
cellular automata. Random Structures & Algorithms 29 (2006) 399-415 

19. Buvel, R., Ingerson, T.: Structure in asynchronous cellular automata. Physica D 1 (1984) 59-68 

20. Huberman, B.A., Glance, N.: Evolutionary games and computer simulations. Proceedings of 
the National Academy of Sciences, USA 90 (Aug. 1993) 7716-7718 

21. Bersini, H., Detours, V.: Asynchrony induces stability in cellular automata based models. In: 
Proceedings of Artificial Life IV, Cambridge, MIT Press (1994) 382-387 

22. Kanada, Y.: Asynchronous ID cellular automata and the effects of fluctuation and randomness. 
In: Proceedings of the fourth conference on artificial life (A-Life IV), MIT Press (1994) 

23. Schonfisch, B., de Roos, A.: Synchronous and asynchronous updating in cellular automata. 
Biosystems 51(3) (1999) 123-143 

24. Fates, N., Morvan, M., Schabanel, N., Thierry, E.: Fully asynchronous behavior of double- 
quiescent elementary cellular automata. Theoretical Computer Science 362 (2006) 1-16 

25. Fates, N., Regnault, D., Schabanel, N., Thierry, E.: Asynchronous behavior of double-quiescent 
elementary cellular automata. In Correa, J.R., Hevia, A., Kiwi, M.A., eds.: LATIN Proceedings. 
Volume 3887 of Lecture Notes in Computer Science., Springer (2006) 455-466 

26. Fuks, H.: Probabilistic cellular automata with conserved quantities. Nonlinearity 17(1) (2004) 
159-173 

27. Fates, N., Morvan, M.: Perturbing the topology of the game of life increases its robustness to 
asynchrony. In Sloot, P.M. A., Chopard, B., Hoekstra, A.G., eds.: 6th International Conference 
on Cellular Automata for Research and Industry Proceedings (ACRI). Volume 3305 of LNCS., 
Springer (Oct. 2004) 111-120 

28. Fates, N., Morvan, M.: An experimental study of robustness to asynchronism for elementary 
cellular automata. Complex Systems 16(1) (2005) 1-27 

29. Rouquier, J.B., Morvan, M.: Combined effect of topology and synchronism perturbation on 
cellular automata: Preliminary results. In Umeo, H., Morishita, S., Nishinari, K., eds.: ACRI. 
Volume 5191 of LNCS., Springer Berlin (2008) 

30. Ermentrout, G.B., Edelstein-Keshet, L.: Cellular automata approaches to biological modelling. 
Journal of Theoretical Biology 160 (1993) 97-133 

31. Silva, H.S., Martins, M.L.: A cellular automata model for cell differentiation. Physica A: 
Statistical Mechanics and its Applications 322 (2003) 555-566 

32. Liggett, T.M.: Interacting Particle Systems, an Introduction. Springer Verlag (2004) 

33. Goles, E., Martinez, S.: Neural and automata networks, dynamical behavior and applications. 
Volume 58 of Maths and Applications. Kluwer Academic Publishers (1990) 

34. Istrail, S.: Statistical mechanics, three-dimensionality and NP-completeness: I. universality 
of intractability for the partition function of the ising model across non-planar lattices. In: 
Proceedings of STOC, ACM Press (2000) 87-96 

35. Rojas, R.: Neural networks: a systematic introduction. Springer (1996) Chap. 13 - The Hopfield 
Model. 

36. Rouquier, J.B., Morvan, M.: Coalescing cellular automata: Synchronization by common random 
source for asynchronous updating. Journal of Cellular Automata 4 (2009) 55-78 



37. Regnault, D.: Directed percolation arising in stochastic cellular automata anaslysis. In: Proceed- 
ings of the 33rd International Symposium on Mathematical Foundations of Computer Science. 
Volume 5162 of LNCS. (2008) 563-574 

38. Lindvall, T.: Lectures on the coupling method. Wiley (1992) 

39. Grimmett, G., Stirzaker, D.: Probability and random process. 3rd edn. Oxford University Press 
(2001) 

40. Grigoriev, S., Priezzhev, V.: Random walk of annihilating particles on the ring. Theoretical 
and Mathematical Physics 146(3) (2006) 411-420 

41. Regnault, D., Rouquier, J.B., Thierry, E.: Stochastic minority on graphs. Technical report, 
LIP, ENS Lyon, Universite de Lyon (April 2008) 

42. Levin, D.A., Peres, Y., Wilmer, E.L.: Markov chains and mixing times. AMS (2009) 

A Characterization of stable configurations on trees 

Let us now study more precisely the structure of the limit set. This characterization allows 
to count the number of attractors on a given tree. 

We first pick an arbitrary vertex of degree 1 and set it as root r in the tree (this introduces 
"parent" and "children" relations). The term "degree" refers to the original graph: a vertex 
of degree three has two children. 

We assign a label to each vertex to count the number of attractors and the size of the 
limit set. The number of attractors will be shown to be the number of acceptable labelings. 
The set of labels we use is {(□,()), (□, 1), (>, 0), ([>, 1)}. V intuitively means "oscillating" 
while "□" means "fixed" (like the recorder symbols play /stop). The second component is 
called the "preferred" state of the vertex. 

Definition 35 (acceptable labeling). A labeling is acceptable if and only if, for each 
vertex i, if the vertex has label 

1. (D,a) then it has strictly more than deg(z)/2 neighbors with label (D,a); 

2. (>,a) then 

(a) if the parent has a label of the form (□,/?), then a = and i has one more child 
labeled (•, 1 — j3) than children labeled (-,/?); 

(b) otherwise, i.e. the parent has a label of the form (>,•), i has one more child labeled 
(•,a) than children labeled (•, 1 — a). 

Note that only vertices of even degree can have a label of the form (>,•). The apparent 
asymmetry in case 2a (imposing a = 0) is there only to avoid double counts in Theorem 37, 
one could as well have defined acceptable labelings with a = 1. 

Theorem 37 shows that a labeling corresponds to an attractor, and Theorem 38 details 
the meaning of a labeling, thus the structure of an attractor. 

For a labeling L, we define a few special configurations: snd(L) is the projection of the 
second component: Li = (-,a) <=> (snd(L))^ = a, and paint (L, a) sets all vertices labeled 



For a configuration c in the limit set, we denote attr(c) the attractor containing c. 



(>, •) to a: 




Lemma 36. If L is an acceptable labeling, snd(L) and paint(L,a) are in the limit set. 

Proof. Consider the configuration paint (L, 0). We consider the vertices in a bottom-up order 
(a vertex is considered after all its children), and update each vertex which is not in its 
preferred state. Updating a vertex makes it goes to its preferred state thanks to the definition 
of acceptable labelings. The current configuration is snd(L). 

Then, we consider the vertices in a top-down order and update those that are in state 
0. When a vertex is considered, its children are in their preferred state, and if the parent 
of a vertex is labeled (>, •) then this parent is in state 1. Again, thanks to the definition 
of acceptable labelings, updating a vertex makes it go to state 1. We get the configuration 
paint (L, 1). 

By symmetry, there is a sequence of updates from paint(L, 1) to paint(L,0). Moreover, 
there is no sequence of updates leading from one of this configurations to change of state of 
a vertex labeled (□, a): the first change of state of such a vertex i would contradict the fact 
that it has more than deg(i)/2 neighbors with label (D,a) and thus in the state a. Which 
implies that, whatever the configuration reached from paint(L,0), there is by monotonicity 
a sequence of updates yielding the configuration paint(L, 1). The full cycle is in the limit 
set. □ 

Theorem 37. Given a tree, there is a bisection between attractors and acceptable labelings. 

Proof. We define a function / that maps an attractor A to an acceptable labeling L. We then 
shows that for every attractor A, attr o snd(/(A)) = A and for every acceptable labeling L, 
/(attr o snd(L)) = L. This imply that / is a bijection. 

To define /, let A be an attractor, we construct L = f(A). For all vertices that have 
the same state a in any configuration of A, define Li := (D,a). All leaves are now labeled. 
We define the labeling of the remaining vertices inductively (in a bottom-up order). Each 
remaining vertex is oscillating, thus has an even degree and an odd number of children. 
Considering a vertex i having all its children labeled: 

— if the parent is already labeled (thus with a label of the form (□,•)), define Li := (>, 0); 

— otherwise, let a be the majority state among the preferred states of the children and 
define Li := (>, a). (Since i has an odd number of children, the majority state is well 
defined) . 



Let us show that the labeling we have just defined is acceptable. With the same argument 
as the proof of proposition 17, the configuration where all oscillating vertices are in state a 
is in A. This configuration is paint(L,a). 

1. Consider a vertex i labeled (□,()) and the configuration paint(L, 1). All the neighbors 
of i not labeled (□,()) are in the state 1. So, there are necessarily more than deg(i)/2 
neighbors labeled (□,()) (and thus in the state 0), otherwise, updating i would make it 
change its state, contradicting the definition of Li = (□, •). By symmetry, point 1 of the 
definition of acceptable labelings is fulfilled. 



2. We show point 2 inductively in a bottom-up order. Consider a vertex i labeled (>, a). 
Let us first show that there is a configuration c G A where all the vertices of the subtree 
having i as root are in their preferred state, except i (which is in state 1 — a). Indeed, 
consider paint(L, 1 — a) and update the vertices not in their preferred state in a bottom- 
up order. Thanks to the induction hypothesis, the label of those vertices is acceptable, 
thus updating them makes them change to their preferred state. This is a sequence of 
updates from a configuration in A leading to c, so c belongs to A. 
We can now show that the label of i is acceptable: 

(a) If the parent has a label of the form (□, •) then Li has been defined as (>, 0). More- 
over, in c, the parent of i is in state j3. In this configuration, updating i cannot make 
the energy decrease, so i must have one more child labeled (-,1 — j3) than children 
labeled (-,/?). So, point 2a is fulfilled. 

(b) Otherwise, the parent is labeled (-,/?). From the definition of L, in configuration c, 
i has a majority of children in state a. Updating it thus makes it change its state. 
So i must have as many neighbors in state a than in state 1 — a. Which means 
that j3 = 1 — a and i has exactly one more child labeled (•, a) than children labeled 
(•, 1 - a). Point 2b is fulfilled. 

Let us now show that for any acceptable labeling V , /(attr o snd(L')) = V . Since 
paint(L / ,0) and paint(L / , 1) are in the attractor attrosnd(L') (cf. Lemma 36), the vertices 
labeled (>, •) in L' are oscillating in this attractor. From the definition of /, these vertices 
are labeled (>, •) in /(attr o snd(L')). 

Recall that there is no sequence of updates leading from snd(L') to change the state of 
a vertex labeled (□, a): the first change of state of such a vertex i would contradict the fact 
that it has more than deg(i)/2 neighbors with label (□,«) and thus in the state a. This 
allows us to conclude that the labelings V and /(attr o snd(L')) have the same cells labeled 
(□,0) and (□, 1), thus they are equal. (Indeed, in the definition of /, the value of a for 
vertices labeled (>, a) is entirely determined by the labeling of vertices labeled (□, •).) 

Finally, let A' be an attractor, we show that attr o snd(/(A')) = A' . f(A') is an acceptable 
labeling so (Lemma 36), paint (/(A'), 1) is in the attractor attr o snd(/(A / )). Moreover, we 
have already noted that the configuration where all the oscillating vertices of A' are in state 1 
belongs to A' (same argument as the proof of Proposition 17). From the definition pf /, this 
configuration is paint(/(A / ), 1). The attractors attr o snd(/(A / )) and A' have the element 
paint (f(A f ), 1) in common, so they are the same attractor. 

This concludes the proof by implying that / is a bijection. □ 

Theorem 38 (Structure of an attractor). Let L be an acceptable labeling. Then for 
every configuration c reachable by a sequence of updates from snd(L), for every vertex i: 

1. If Li = (D,a) then q = a (this is why intuitively means "fixed"). 

2. If Li — (>, a) (in this case deg(i) is even: i is not the root, which has degree 1) then 

(a) if the parent has a label of the form (□, i is in the state appearing in majority 
among its neighbors (no constraint in case of equality); 

(b) otherwise 



— ifiis in its preferred state a, its children labeled (•, a) are in their preferred state 
a 

— otherwise, all its children not in their preferred state are in the same state as i 
(the state 1 — a). 

Proof Recall that there is no sequence of updates from snd(L) to the firing of a vertex 
labeled (D,a). This ensures point 1. 

Since snd(L), thus c, are in the limit set, the energy cannot decrease. It follows that any 
vertex i has always at least deg(z)/2 neighbors in the same state as i. This proves point 2a. 

We prove the last point (2b) by recursion. Note that the assertion of the theorem is 
clearly true for snd(L). Let c be a configuration verifying the assertion, c' a configuration 
reached from c by updating a vertex z, it is sufficient to prove that the assertion is true for 
c' . Precisely, we prove that it is true for i and its neighbors. 
For i: 

— If i is in its preferred state a in c, then its children labeled (there are deg(z)/2 
such children thanks to the definition of acceptable labelings) are then in state a in c 
and thus in d . Which means that i can change its state only if all other children are in 
state 1 — a. Point 2b stays true for i in d . 

— Otherwise, it is either inactive and point 2b stays true in c', or active. In the latter 
case, i has as many neighbors in each state (because c is in the limit state), all children 
of i not in their preferred state are in the same state as i (by recursion hypothesis) 
and i has one more child labeled than children labeled (-,1 — a) (definition of 
acceptable labelings). These three conditions imply that all children of i are actually in 
their preferred state. So, d fulfills point 2b. 

For a neighbor w of i: 

— If w is a child of z, nothing has changed for the sons of w: point 2b stays true. 

— If w is the parent of i, it is sufficient to consider the case where i is in the same state a 
as w (there is no condition to check in the other case). 

• If a is the preferred state of i then i has at least deg(z)/2 + 1 neighbors in state a: 
its children labeled and its parent w. So i is inactive. 

• Otherwise, i and w are not in their preferred state, and point 2b stays true, whether 
or not i change its state. □ 

B Omitted proofs: Bound on the Hitting Time of a 2D Finite 
Markov Chain 

B.l Background on Markov chain theory 

We recall only the necessary background on Markov chains to get a bound on the hitting 
time of a 2D finite Markov chain. For a gentle introduction and proofs, we refer for instance 
to Chapter 10 of [42]. 

Let (X t )teN be a Markov chain. We note r& the hitting time of 6, i.e. the first time the 
Markov chain is in state b: 

Tb := min{t ^ | X t = b} . 



If (X, P) is a Markov chain reversible with respect to the probability 7r, the conductance 
of an unoriented edge (x, y) is 

g(x, y) := tt(x)P(x, y) = n(y)P(y, x) . 

Let a and b be two distinguished vertices, representing source and sink. We use the following 
potential, or voltage, of a vertex: 

V{x) :=F(r a <n\X = x) . 

Clearly V(a) = 1 and V(b) = 0. Now, define the current flow on oriented edges as 

I(x, y) := g(x, y) (V(x) - V(y)) and ||/|| := £ I(a, x) . 

X 

The effective resistance between a and b is 

Theorem 39 (Commute time identity). 

E (r b | X = a) + E (r a | X = 6) = ^(a, 6) 

In our case, g = 1, so 7£(a, 6) is an upper bound for the average hitting time E (t& | Xo = a). 
Here is how one can bound 7£(a, b). A /tow from a to b is a function on oriented edges which 
is antisymmetric: 0(x,y) = —6(y,x) and which obeys Kirchhoff's vertex law: 

Vv £ {a,b} ^0(x,v)=0 . 

X 

This is just the requirement "flow in equals flow out". The strength of a flow is 

||0|| 

X 

Theorem 40 (Thomson's Principle). For any finite connected graph, 
1Z(a, b) — inf {E{9) \ a unit flow from a to b} 

where E(6) := V (g( f '^f . 



B.2 Application to our case 

If a neighbor does not exist (because the vertex is on the border), the edge points to the 
vertex itself. Each edge has the same probability \. 



An invariant probability is the uniform probability ir : x H> Our Markov chain is 
reversible with respect to this probability. So we can use the definitions of section B.l: 



\/x,y g(x,y) 



4nm 



Thanks to Thomson's principle, it is sufficient to construct a flow from a to b to get an 
upper bound on 7Z(a,b). If a = (i,j) and b = (n, m), we consider the trivial (and far from 
optimal) flow 

'9((kJ), (k + 1, j)) := 1 if k < n 
(9((n, fe), (n, k + 1)) := 1 if j < fc < m 
^(x,2/):=0 elsewhere. 

That is, a flow on a single path from a to b. 

E(0) ^ (n + m)4nm. We conclude with the commute time identity that the average 
hitting time is 0(n 3 ), assuming w.l.o.g. that n ^ m. 



